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The gluon propagator is computed in large asymmetric lattices, accessing momenta around 100 
MeV and smaller. Our study tries to check the compatibility of the recent solutions of the gluon- 
ghost Dyson-Schwinger equations with lattice results. In particular the exponent K, which char- 
acterizes the solutions for the infrared, is measured. Results favours a vanishing zero momentum 
gluon propagator. We also report on the compatibility of various functional forms used to fit the 
Dyson-Schwinger solution for the full range of momenta and the lattice gluon propagator. 
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Quantum chromodynamics (QCD) is the theory that describes the interaction between quarks 
and gluons. Two first principles approaches to non-perturbative QCD are Dyson-Schwinger equa- 
tions (DSE) and the lattice formulation of QCD (LQCD). Solving the DSE, an infinite tower of 
equations, requires a truncation scheme and a parametrization of higher order vertices [p. On 
the other hand, in LQCD one has to care about finite volume [||] and finite space effects. If han- 
dling finite space effects requires smaller lattice spacing, i.e. higher /3 values, finite volume effects 
requires large volumes. 

Recently, there has been some progress on the solution of the DSE. Assuming ghost domi- 
nance, the equations were solved in the infrared (IR) limit. The solution shows a vanishing zero 
momentum gluon propagator [Q]. Moreover, the coupled gluon-ghost equations were solved nu- 
merically and analitical functions found that fitted the solution [|], Q]. On the lattice, the gluon 
propagator has been revisited a number of times. However, due to the lattice sizes used in the 
investigations, the IR region was not yet properly accessed. The motivation of our study is to look 
at the deep IR gluon propagator and to try to check the compatibility between the two approaches, 
i.e. to test the IR DSE solution for the gluon propagator and to see if the same functional forms 
describe the lattice data and the DSE solution. In order to try to investigate the deep IR region, we 
consider large assymetric 4D lattices, i.e. lattices with a very large time extension [||]. 

1. The lattice setup 

In this investigation we consider SU(3) pure gauge, Wilson action at j3 = 6.0. The value of 
the lattice spacing to convert lattice into physical units being a~ l = 1.943 GeV (J7|. For the lattices 
16 3 x 128 and 16 3 x 256, the gauge configurations 1 were generated as described in For the 
smallest (largest) lattice 164 (155) configurations were generated. The large time extension allows 
to access momenta of about 95 MeV for the smallest lattice and 48 MeV for the largest lattice. For 
the definition of the gluon propagator, notation and the gauge fixing method see [gj. 

2. Gluon dressing function 

In figure [T] the bare gluon dressing function is shown for all pure temporal and pure spatial 
momenta. The data shows volume effects, with the spatial momenta giving larger values for the 
dressing function. The data suggests that one should not mix temporal and spatial momenta. 

For the IR study, we will consider only pure temporal momenta. For these momenta and 
the two lattices, the gluon dressing function and the propagator are reported in figure |2[ The two 
data sets are compatible within errors. Note that this does not mean that finite volume effects are 
negligable [§]• 

2.1 IR gluon dressing function 

The IR analytical solution of the DSE for the gluon dressing function being Z(q 2 ) ~ (q 2 ) 2K , 
we first try to fit the lattice data to this expression 2 . An acceptable % 2 jd.o.f . = 0.40 was obtained 

'All configurations were generated with MILC code http : / /physics .indiana.edu/~sg/milc. html. 
2 The fits were done separately for each lattice. The statistical errors on all fitting parameters were estimated using 
1600 bootstrap samples. 
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Figure 1: Gluon dressing function for the smaller (left) and larger (right) lattices. 
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Figure 2: Gluon dressing function (left) and propagator (right) for pure temporal momenta. 



only by fitting the first three lowest momenta (excluding the zero momentum) for the largest lattice 
(\q\ < 143 MeV), meaning that the DSE solution seems to be valid only for momenta below 150 
MeV. The measured K = 0.4859^24 does not support a vanishing zero momentum gluon propagator. 
In order to test the validity of the result, corrections to the DSE solution were considered, fitting 
the smallest possible range of momenta. The results are in table 1. 



16 3 x 256 


9max 


K 




X 2 /d.o.f. 


(0(q 2 ) 2K (l+aq 2 ) 


\9\MeV 


0.5070: 


-36 
-50 


0.44 


(0(q 2 ) 2K {\+aq 2 + bq 4 ) 


238MeV 


0.5131: 


-67 
-64 


1.03 


(0(q 2 ) 2K {l+aq 2 + bq 4 + cq 6 ) 


286MeV 


0.5146: 


-70 
-96 


1.36 



Table 1: Fits to polinomial corrections to the DSE analytical solution. 



It is interesting to observe that, with corrections to the analitycal DSE solution, K becomes 
larger than 0.5. Moreover, the results reported in table 1 are compatible within errors and give an 
error weighted mean value of K = 0.5108(13). Note that the K just reported are not compatible 
with the direct fit to the DSE solution. 
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In [Q] the DSE solution was well described by the two expressions 



z cu t(q 2 



ft) 



q 



A QCD + Q 2 



-'pole { 



ft)- 



2k 



QCD 



2 k 



(2.1) 



+ {q 2 



which associate, respectively, a cut and a pole to the propagator in the IR region. The lattice data is 
well described by both expressions with no clear preference for the cut or the pole. The results of 
the fits, for the momenta region starting at the lowest nonzero momentum up to q max , are reported 
in table |2[ Note that all fits have good % 2 jd.o.f . and, for each function, the fitted parameters are 
compatible within one standard deviation. 





Lattice 


qmax 


K 


Aqcd (MeV) 


X 2 /d.o.f. 


^cut 


16 3 x 256 
16 3 x 128 


664MeV 
510MeV 


0.5090+^ 
0.51171$ 


409+4 
417±| 


0.71 
1.25 




16 3 x 256 
16 3 x 128 


664MeV 
510MeV 


0.5077+^ 
0.5100+j? 


409+* 
416+j 


0.69 
1.15 






Table 2: 


Z cut and Zp i e 


fits. 





2.2 UV gluon dressing function 

As in previous studies [0, |8]] we fit the UV data to the one-loop perturbative gluon dressing 
function, 




(2.2) 



where y = —13/22 is the gluon anomalous dimension. In the UV region two types of momenta 
were considered, namely pure temporal and cylinder plus conic cuts as defined in [§]. The results 
can be seen in table || for the largest acceptable momentum range. They show that the UV region 
is well described by the one-loop perturbative result, with the temporal momenta following the 
perturbative behaviour earlier than the cylindrical+conical cut momenta. Note that since the fits are 
to the bare data, the fitted A's should not be compared directly with standard values. 



Lattice 


Temporal 

l.SGeV < \\q\\ <3.9GeV 
A X 2 /d.o.f. 


Cuts 

l.SGeV < \\q\\ <7.SGeV 
A X 2 /d.o.f. 


16 3 x 128 
16 3 x 256 


1006+^MeV 1.2 
990+jMeV 1.1 


640+jMeV 1.1 
642>t\MeV 1.2 



Table 3: UV fits. 



2.3 The gluon dressing function 

In what concerns the full data set, it was tested against the following formulas Zfi t (q 2 ) = 
B{q 2 ) ((Xfit(q 2 )) 7 , with B(q 2 ) = Z po j e ^ cut {q 2 ) and using two different definitions for the running 
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coupling, 

a P {q 2 ) 

<W<7 2 ) 



r 

<P — 



o(0) + 



471 



(. 



A QCD A> ^H^/^QCd) 



<1 2 / A2 QCD ~ 1 



a(0) 



In 



e+fli^/Alo,)* 



,A) = H, (2-3) 



(2.4) 



proposed in [Q |5p. The fits with (Xp{q 2 ) had always % 2 jd.o.f . > 2, i.e. the lattice data is not 
described by such a running coupling constant. The lattice data adjusts better to CCtNiq 2 ) as can be 
seen by the results reported in table 4. Note that the K values are now larger and closer to the DSE 
figures, than those measured in the IR region. 



B = 




K 


A QCD (MeV) 


a\ 


a 2 


X 2 /d.o.f. 


16 3 
16 3 


x 128 
x256 


0.5435+j^ 
0.5244+^ 


364+4 
374^2 


0.0062+ 3 
0.00721 3 


2.44+j 
2.42411° 


1.82 
1.73 


B = 


Z pole 


K 


K QCD (MeV) 


a\ 


ai 


X 2 /d.o.f. 


16 3 
16 3 


x 128 
x256 


0.5335+J 3 
0.5217+J4 


373+| 
377+1 


0.0081 I 3 
0.0082+} 


2.36+J 
2.36+} 


1.65 
1.61 



Table 4: Fits to all lattice data. 

The value of the zero momentum running coupling can be computed from the asymptotic 
behaviour of the QCD /3-function and from (Xln ip ) > oc{0) = (47r//3o)c?2- Note that the values 
reported in table 5 are not too far from the DSE numbers. 



o(0) 






16 3 x 128 


2.791 2 


2.70+J 


16 3 x 256 


2.771} 


2.701} 



Table 5: Values of a(0) extracted from (Xln(p 2 ). 



3. Gribov copies 

For the smallest lattice, we look for possible Gribov copies effects comparing the method used 
in [^] with that described in JTo|]. The ratio of the propagators computed from the two methods, 
Djd/Dceasd, for pure temporal momenta is shown in figure ||. Although we have a ratio that, 
for many momenta, is not compatible with one, we cannot conclude in favour of any systematic 
effect 3 . This means either that our statistics is too small to resolve Gribov copies or they do not 
play a significant role for such lattices. 

3 In M the effect of the Gribov copies was estimated as a 2 to 3 c effect. 
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Figure 3: Gribov copy effects in 16 3 x 128 configurations. 



4. Conclusions 

The lattice data seems to be compatible with functional forms that describe the gluon-ghost 
DSE numerical solution. In what concerns the vanishing of the zero momentum gluon propagator, 
our data does not provide a clear answer but it favours a vanishing propagator in agreement with the 



DSE solution, flow equations studies [11, 12] and time independent stochastic quantisation [13]. 
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